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Abstract 

Using exact diagonalization method,we studied the ground state of the 
anisotropic molecular magnets and find a critical point in the transverse 
field,which may divide the quantum tunneling region into two different parts. 
Possible ways to observe and take advantage of this point by varying the 
transverse field are suggested. 
Recently the macroscopic quantum phenomena in large spin anisotropic magnets is 
attracting much attention among both theorists and experimentalists [1-2]. Many efforts 
has been done to measure the step-wise magnetization loop, classify different temperature 
regimes, and find different theoretical explanation for each. Quantum resonant tunneling is 
the underlying mechanism that leads to abrupt change in the magnetsation. However, it's not 
easy to take into account various factors that may influence the result quantitatively, which 
sometimes results in sheerly different pictures [3-6]. Fortunately,the model Hamiltonian has 
been established and widely used in theoretical analysis. Many work has been starting from 
the single large spin model, treating other parts as mean fields,with some static or dynamic 
distribution function. There are three different temperature regimes that can be described as 
different phases:thermal-activated regime, thermal-assisted regime, and purely quantum tun- 
neling regime. Apparent phase transition(PT) has been studied by analytical and numerical 
means and both first and second order phase transition are found [7]. In this paper, thro ugh 
a properly denned order parameter, we find a kind of quantum phase transition(QPT) for 
Mni 2 which is induced by continuous change of transverse field in the third regime. Its 
experimental consequences is suggested. 
Firstlythe model Hamiltonian of Mni2 is 

H = -D(S Z ) 2 - B(S Z ) 4 - C((S + ) 4 + (S-) 4 ) - gfi B h x S x - gfi B h z S z (1) 

where the first and second terms are diagonal part of the crystal field,the third term 
is the off-diagonal part. The last two terms represent external field. In writing down it, 
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we neglect all the other internal degrees of freedom of the cluster on one lattice site,i.e.,the 
high energy spectrum which can be deviating the description of fixed spin magnitude, for 
the antiferromagnetic coupling between internal spins (4 spins 3/2 and 8 spins 2 for Mni 2 ) 
is rather large. This is called the "mesoscopic " approach. The parameters 

are determined by various measurements. In the following we choose D m 0.60K, B/D 
~ 0.002, C ~ 3 • 10~ 5 fT, which are obtained in EPR experiments. Furthermore we set h z = 
to study pure transverse field effect. 

In the pure quantum relaxation regime,without the longitudinal field,tunneling occurs 
between degenerate ground states —S and +S. In the Landau-Zener theory the transition 
probability is given by: 

P_ 5i5 = 1 - exp(-r_ Sj5 ) (2) 

with 

= irA 2 { _ SiS) /[4hS gf jL B dH/dt] (3) 

which is determined by tunneling splitting as well as sweeping rate of applying field. 
In this paper, however, we will study a related problem — the characteristic of the ground 
state — and understand the physics in a different way. We consider the states with largest 
| (S z ) |, denoted as \g) , which is also energetically lowest in most cases. Generally \g) can be 
expanded as 

s 

\g)= Yl C ™W (4) 

m=-S 

Now we define order parameter as the difference of probability of the \±S) components 
in \g) : 

5 = \c s \ - \c- s \ (5) 

For situations where quantum tunneling due to ground state dynamics dominates in 
sufficiently low temperature, 8 measures the ability of ground state \g) to transfer spins 

3 



from one side in magnetic spectrum to other side. Some cases corresponding to (1)B = 

0, C = 0.(2)5 = 0.002, C = 0.(3)5 = 0,C= 10" 5 ,are studied and demonstrated in Figure 

1. The 5 line changes abruptly near some point of h x , showing a quantum critical point h c 
in the classical limit S — > oo. Across this narrow crossover region,the difference of \±S) 
components in \g) (i.e., S), changes from almost 1 to O.When \h x \ > h c , numerical results show 
equal components of the "bare" states S = ±10, in our case. The effect of B and C term 
on the 5 line,especially on the crossover position are easy to understand. B term tends to 
strengthen the anisotropyleading to higher barrier between the positive and negative half of 
the magnetic spectrum. Thus a larger value of \h x \ is needed to cause tunneling combination. 
However, the C term is transverse and will help tunneling, and because of its four-ordered 
operator form,can influence the results greatly by a little change of its value(10 -5 relative to 
D, and 10~ 2 to B). The combined effect is shown in Figure 2. The C term quickly decreases 
the difference of the two "phases". In some systems, C ~ 3 • 10~ 5 so maybe we need some 
trick to really observe the critical point. 

Also we calculated the distribution on bare states of \g).As is evident from Figure 3. The 
|±S") components are sizable near the crossover region in the 5 = phase. The experimental 
consequence is, in a well prepared sample, if we drive all sample sites to \S = 10) state, a 
considerable portion of lattice site can tunnel to the \S = —10) state through the ground 
state channel, if the transverse field is properly controlled. Or, by increasing the transverse 
field across the crossover region,we can "pump" many sites from one side of the magnetic 
spectrum to another. Such scenario is related to relaxation experiment ||. 

To conclude,by means of exact diagonalization,we find a crossover in the Hamiltonian of 
molecular magnets,through which the difference of amplitudes of \±S = 10) states in ground 
state are changed from almost 1 to zero. This property might have some experimental 
consequences in the purely quantum regime of the magnetization relaxation which is caused 
by tunneling through the ground state channel. 
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Caption: 

Figure 1. The amplitude difference 5 of the \±S) components in the ground state \g) 
for three groups of B, C values. The 5 line changes abruptly near some h x points. The effect 
of B is to broaden the cental "phase" of the diagram, while the C value effects the opposite 
way. 

Figure 2. For larger C values the central phase narrows quickly and will disappear at 
large enough C value. This made observation of the critical point of transverse field more 
subtle and difficult. Here we plotted two lines corresponding to C — 1CT 5 and C — 2 • 1CT 5 . 

Figure 3. The amplitudes of the ground state projected to bare states:]— S) \S), for 
cases (a) h x — 1, (6) h x = 2.3 ,(c) h x = 2.7, (d) h x = 3.3. Both B and C are choosen to be 
zero. The distribution tends to become symmetric as h x transcends some value. 
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